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Two-dimensional  arrays  of  coupled  nanomechanical  plate-type  resonators  were  fabricated  in  single 
crystal  silicon  using  e-beam  lithography.  Collective  modes  were  studied  using  a  double  laser  setup 
with  independent  positioning  of  the  point  laser  drive  and  interferometric  motion  detector.  The 
formation  of  a  wide  acoustic  band  has  been  demonstrated.  Localization  due  to  disorder  (mistune) 
was  identified  as  a  parameter  that  limits  the  propagation  of  the  elastic  waves.  We  show  that  all  400 
resonators  in  our  20  X  20  array  participate  in  the  extended  modes  and  estimate  group  velocity  and 
density  of  states.  Applications  utilizing  the  resonator  arrays  for  radio  frequency  signal  processing  are 
discussed.  ©  2006  American  Institute  of  Physics.  [DOI:  10.1063/1.2190448] 


The  collective  behavior  of  coupled  nanoelectromechani¬ 
cal  system  (NEMS)  resonators  can  largely  enhance  perfor¬ 
mance  of  resonator-based  devices  [radio  frequency  (RF) 
filters,1  added  mass  sensors,2  magnetometers,3  etc.]  while 
preserving  the  most  attractive  features  of  microelectrome¬ 
chanical  resonators — superior  quality  factor,  microminiature 
size,  and  compatibility  with  complementary  metal-oxide 
semiconductor  (CMOS)  processes.  When  the  size  of  the  ar¬ 
ray  exceeds  hundreds  or  thousands  of  resonators,  new  capa¬ 
bilities  can  emerge  that  are  unavailable  from  a  single  NEMS 
device.  A  coupling  between  neighboring  resonators  creates 
an  artificial  crystal — a  manmade  analog  of  a  crystal  lattice 
described  in  solid  state  physics  in  terms  of  the  tight  binding 
approximation.4  An  acoustic  band,  formed  by  splitting  indi¬ 
vidual  resonant  modes,  opens  the  opportunity  for  operating 
with  broadband  excitations  implemented  as  elastic  waves 
propagating  through  the  array.  The  theoretical  language  de¬ 
veloped  for  closely  related  electromagnetic  systems:  coupled 
microwave  cavities,5  traveling  wave  amplifiers,6  and  optical 
slow  wave  structures7  (SWS),  can  be  largely  adopted  to  de¬ 
scribe  the  elastic  waves  in  mechanical  resonators  arrays.  A 
rich  behavior  with  a  collection  of  phenomena  ranging  from 
intrinsic  localized  modes8  (ILM)  to  solitons9  can  be  exhib¬ 
ited  by  such  systems. 

In  practice,  however,  fabrication  of  these  arrays  appears 
to  be  challenging  due  to  micromachining  tolerances  that 
cause  random  mistuning  of  individual  NEMS  resonators.  For 
larger  structures  (270  yum  long  clamped-clamped  beams, 
/reS  =179  kHz)  the  frequency  spread  can  be  confined  within 
0.8%  and  existence  of  collective  modes  was  reported  for  a 
linear  array  of  67  of  such  resonators.10  However,  the  problem 
becomes  more  severe  when  the  scaled  down  devices  have 
nanometer  dimensions  and  exhibit  RF  frequencies.  Excessive 
mistuning  (referred  to  as  a  disorder)  leads  to  localization  of 
the  acoustic  excitations  in  the  array  in  a  way  similar  to 
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Anderson’s  localization11  (AL)  and  raises  the  question 
whether  extended  modes  can  be  realized  in  large  arrays  of 
coupled  nanomechanical  resonators. 

In  order  to  suppress  the  localizing  effects  caused  by  the 
disorder  we  expanded  the  dimensionality  of  the  arrays  from 
linear  to  two  dimensional.  According  to  AL  theoretical  mod¬ 
els,  in  two-dimensional  (2D)  structures  the  elastic  waves  are 
still  expected  to  be  localized  for  an  infinitesimal  amount  of 
disorder;  however,  the  localization  length,  \,  can  be  ex¬ 
tended  by  orders  of  magnitude  compared  to  one-dimensional 
(ID)  case.13'14  In  this  letter  we  present  our  results  on  the 
fabrication  and  experimental  study  of  collective  elastic  exci¬ 
tations  in  2D  arrays  of  coupled  plate-type  nanomechanical 
resonators.  Each  resonator  in  the  array  is  comprised  of  a 
silicon  square  plate  (2X2  /xnr  area,  45  nm  thickness)  sup¬ 
ported  by  a  square  Si02  pillar  at  the  center.  The  coupling 
between  neighbors  is  provided  by  narrow  beams  (200  nm 
wide,  45  nm  thick)  connecting  the  centers  of  the  square  sides 
(Fig.  1).  Commercially  available  silicon  on  insulator  (SOI) 
wafers  [Si  layer,  45  nm;  buried  oxide  (BOX),  150  nm]  were 
used  in  fabrication.  The  geometry  of  the  arrays  was  defined 
by  e-beam  exposure  of  hydrogen  silsesquioxane  (HSQ)  resist 


FIG.  1.  SEM  image  showing  a  section  of  the  paddle  resonator  array.  The 
size  of  the  silicon  square  paddles  is  2  X  2  yum2,  with  a  thickness  of  45  nm. 
The  Si02  pillar  supporting  the  center  of  the  paddle  has  a  cross  section 
0.65  X  0.65  /mi2  and  is  150  nm  tall.  The  connecting  bars  are  2  yum  long  and 
200  nm  wide. 
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FIG.  2.  A  vibration  spectrum  and  mode  shapes  for  a  single  paddle  resonator 
(size  identical  to  that  shown  in  Fig.  1)  with  connecting  arms  coupled  to  a 
frame. 

(Dow  Corning  FOX-12;  thickness,  100  nm)15  spun  on  the 
top  of  the  SOI  layer.  After  developing  the  resist,  a  dry  etch  in 
a  mixed  C4F8  and  SF6  plasma  with  the  FOX-12  layer  as  a 
mask  was  used  to  transfer  the  pattern  onto  the  Si  layer.  A 
subsequent  wet  etch  in  concentrated  FIF  (48%)  removes  the 
rest  of  the  HSQ  and  partially  undercuts  the  BOX  layer  re¬ 
leasing  the  periphery  of  the  Si  square  plates.  The  wet  etch  is 
timed  in  order  to  leave  a  Si02  supporting  pillar  at  the  center 
of  the  resonator. 

To  study  the  vibratory  response  of  the  array,  we  use  an 
optical  setup  similar  to  one  described  elsewhere.16  A  point- 
type  excitation  is  provided  by  modulation  of  the  intensity  of 
the  blue  diode  laser  (wavelength,  412  nm)  focused  on  a 
single  paddle  resonator  within  the  array.  The  resulting  mo¬ 
tion  can  be  detected  by  focusing  another  laser  (red  cw  HeNe, 
633  nm)  at  any  point  of  the  array  and  measuring  the  modu¬ 
lation  of  the  intensity  of  the  reflected  light.  This  modulation 
is  provided  by  a  Fabry-Perot  interferometer  formed  by  the 
paddle  itself  (acting  as  a  moving  semitransparent  mirror)  and 
silicon  substrate  underneath.17  Low  absorption  in  the  very 
thin  (45  nm)  Si  layer  allows  us  to  increase  the  laser  power  up 
to  5  mW  without  damaging  the  nanostructures.  Independent 
positioning  of  the  blue  (drive)  and  red  (detect)  lasers  pro¬ 
vides  the  flexibility  that  is  necessary  to  study  the  propagation 
of  the  elastic  waves  through  the  array. 

Figure  2  shows  the  resonant  spectrum  of  an  isolated 
single  unit  of  the  array — a  test  structure  comprised  of  a 
paddle  resonator  with  2  /uni -long  connecting  arms  bound  to 
an  undercut  frame.  Following  the  results  of  finite  element 
analysis  (FEA),  we  denote  low  frequency  peaks  (45-55 
MHz)  )  as  bridge-type  vibrations  of  coupling  arms.  A  high 
frequency  group  of  resonances  (80-100  MHz  range)  origi¬ 
nates  from  the  out-of-plane  motion  of  the  paddle  itself  and 
includes  y0o  and  yu  modes.18  All  the  peaks  show  a  quality 
factor  Q  ~  3000  in  vacuum  ~10-7  Torr. 

A  drastic  transformation  of  the  vibration  spectrum  is  ob¬ 
served  when  the  paddle  resonators  are  coupled  in  a  square 
array.  The  inset  in  Fig.  3  shows  that  in  a  20  X  20  array 
(2X2  /zm2  paddles,  2  /urn  long  connecting  arms,  0.67  /zm 
undercut)  the  individual  levels  of  the  resonators  split  and 
create  a  27  MHz  wide  band  (91-118  MHz)  exhibiting  hun¬ 
dreds  of  sharp  peaks.  This  spectrum  was  acquired  with  the 
blue  laser  providing  excitation  at  the  center  of  the  array  and 
the  red  laser  focused  on  a  resonator  20  /zm  (five  cells)  away 
from  the  driving  point. 

In  order  to  study  the  propagation  range  of  the  acoustic 
waves  over  our  array,  spectra  similar  to  that  shown  in  the 

inset  to  Fig.  3  were  acquired  at  various  distances  from  the 
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FIG.  3.  The  inset  shows  an  acoustic  band  in  the  vibration  spectrum  of  the 
20  X  20  array  of  coupled  resonators.  The  excitation  was  provided  at  the 
center  of  the  array  and  the  detection  laser  was  positioned  5  unit  cells  away. 
Solid  circles  in  the  main  plot  demonstrate  spatial  decay  of  the  average 
acoustic  energy  of  vibrations  as  a  function  of  separation  distance  between 
the  drive  and  detect  laser  spots.  The  dotted  line  shows  an  exponential  fit  for 
the  decay. 

driving  point  (excitation  was  always  provided  at  the  center  of 
the  array).  The  solid  line  in  Fig.  3  shows  the  spatial  decay  of 
the  band-averaged  acoustic  energy  Ea,  calculated  as  a 
squared  amplitude  of  the  motion  of  a  particular  resonator 
(separated  by  distance  R  from  the  driving  point)  integrated 
over  91-118  MHz  frequency  band.  The  exponential  decay 
of  Ea(R )  provides  a  crude  estimate  of  X  ~  32  /zm  (eight  lat¬ 
tice  sites)  for  the  localization  length  calculated  according  to 
Ref.  12, 

2  In  t(r,r'-E) 

—  =  -  lim  — : - - — , 

^  V~r'\ — >°°  r  ^  I 

where  t{r,r'  \E)  is  the  energy  transmission  probability.  The 
high  quality  factor  of  our  MEMS  resonators  allows  us  to  rule 
out  dissipation  as  the  factor  that  limits  propagation  of  acous¬ 
tic  waves  in  array.19  We  attribute  the  spatial  decay  of  ampli¬ 
tude  to  disorder-driven  localization,  determined  by 
constructive/destructive  interference  of  the  acoustic  waves 
that  follow  different  paths  along  the  2D  array.  A  theoretical 
value  for  the  localization  length  X  can  be  deduced  from  the 
known  value  of  the  W/V  parameter.11'20  We  estimate  the  dis¬ 
order  W  as  a  doubled  standard  deviation  (averaged  over  dif¬ 
ferent  modes)  of  the  frequency  spread  crstd=3.6  MHz  exhib¬ 
ited  by  our  test  structures  (single  resonators).  The  coupling 
parameter  V  was  extracted  from  the  width  of  the  yu  band 
(A/ 1 ,  =4  MHz  as  a  result  of  FEA  calculations).  The  resulting 
value  W/V=  7  provides  an  estimated  Xtheor~  18  unit  cells  ac¬ 
cording  to  the  results  of  numerical  simulations  for  2D  disor¬ 
dered  media.14  We  find  our  experimental  X  value  to  be  in 
reasonable  agreement  with  Xtheor,  considering  the  limited  size 
of  our  microfabricated  arrays  and  fact  that  we  have  used 
bandwidth  A fu  estimated  for  one  particular  mode,  while  our 
acoustic  band  is  generated  by  at  least  three  overlapping 
modes  shown  in  the  high  frequency  group  in  Fig.  2. 

Stronger  coupling  and  increased  density  of  states  for  the 
collective  modes  of  the  array  can  enhance  the  propagation 
range.  In  our  2D  arrays,  we  can  implement  both  approaches 
by  tailoring  the  geometry  of  the  array.  Shortening  the  con¬ 
necting  arms  provides  obvious  enhancement  of  the  coupling 
strength  between  plate-type  resonators.  Shrinking  the  central 
Si  CL  pillar  supporting  the  square  resonators  (by  prolonging 
HF  wet  etch  time)  reduces  frequencies  for  the  plate-type  vi¬ 
brations.  By  combining  the  effects  of  shortening  the  arms 
to  AIP  license  or  copyright,  see  http://apl.aip.org/apl/copyright.jsp 
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FIG.  4.  Enlarged  portion  of  the  vibration  spectrum  and  phase  shift  plot  for 
20X20  paddle  array  with  stronger  coupling  (paddle  size,  2X2  gnf :  con¬ 
necting  arms,  1  jitm)  and  wider  undercut  (Si02  pillar,  0.54  X  0.54  /xm2)  The 
spectrum  was  acquired  with  a  separation  of  40  /im  between  drive  and  de¬ 
tection  points.  The  slope  of  the  A<pi  to)  line  corresponds  to  a  group  velocity 
Vg, — 10  m/s. 

and  enhancing  the  plates’  undercut,  one  can  shift  the  frequen¬ 
cies  of  the  bridge-type  modes  for  the  connecting  arms 
(45  MHz  cluster  on  Fig.  2)  up  to  the  limit  where  they  start 
overlapping  with  the  plate-type  modes  of  the  paddles,  thus 
increasing  the  density  of  states.  A  rectangular  20  X  20  array 
with  1  /am  long  arms,  2X2  /am2  paddles,  and  0.73  /am  un¬ 
dercut  was  fabricated  according  to  this  recipe.  We  have 
found  that  the  spatial  decay  in  the  amplitude  of  the  vibration 
over  the  entire  20  X  20  array  was  only  of  the  order  of  50%, 
i.e.,  all  400  resonators  appear  to  be  involved  in  the  collective 
motion. 

By  counting  resonance  peaks  in  the  vibration  spectrum 
of  this  array  (part  of  the  spectrum  is  shown  in  Fig.  4),  one 
can  estimate  density  of  states  as  13  peaks  per  MHz  (to  be 
compared  with  —30  MHz-1  deduced  from  FEA  calcula¬ 
tions).  The  unwrapped  phase  dependence  shown  in  Fig.  4 
provides  a  mean  value  for  group  delay, 

Tdelay  =  dip/ d(x)  —  5.8  faS  . 

From  the  dimension  of  the  array,  L= 60  /am  one  can  estimate 
group  velocity,  Vgr~10  m/s.  Pronounced  kinks  in  the 
A  <p(co)  graph  at  67  and  71  MHz  illustrate  that  the  group  ve¬ 
locity  can  have  a  strong  frequency  dependence.  In  general, 
chromatic  aberration  can  alter  the  group  velocity  by  orders  of 
magnitude  in  arrays  of  coupled  resonators21  and  may  be  ex¬ 
ploited  for  RF  signal  processing.  By  converting  RF  pulses 
into  elastic  waves  propagating  through  an  inhomogeneous 
dispersive  array,  one  may  be  able  to  implement  for  example, 
a  broadband  spectrum  analyzer  mimicking  Newton’s  optical 
prism.  More  elaborate  methods  might  exploit  both  normal 
and  anomalous  dispersions  and  effects  based  on  nonlinearity 
(wave  mixing,  soliton  formation,  etc.). 

As  the  next  step,  we  are  planning  to  study  diffusivity  and 
dispersion  in  our  structures  by  fabricating  larger  arrays  and 
measuring,  in  the  time  domain,  propagation  of  short  wave 
packets.  Currently  we  are  also  taking  measures  to  further 
refine  the  fabrication  process  and  reduce  the  disorder. 

To  conclude  we  have  fabricated  and  tested  2D  arrays  of 
coupled  nanomechanical  resonators  and  demonstrated  the 


existence  of  extended  modes.  Localization  of  the  elastic 
waves  due  to  fabrication-induced  disorder  is  a  key  technical 
challenge  but  it  can  be  addressed  by  tailoring  the  array  de¬ 
sign.  We  demonstrate  that  all  400  resonators  in  our  array 
were  involved  in  the  acoustic  band  formation  and  give  an 
estimate  of  13  MHz-1  for  the  density  of  states.  We  anticipate 
that  NEMS  arrays  may  ultimately  have  a  wide  range  of  ap¬ 
plications  in  sensing  and  RF  signal  processing. 
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